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SUMMARY 


This  report  describes  a  continuation  of  the  work  reported  in  •»  ■  /I’ I Af 
Technical  Note  by  A.  C.  Campbell  and  B.  Etkin  titled  "The  Response  of  Cylindri¬ 
cal  Structure  to  a  Turbulent  Flow  Field  at  Subcritical  Reynolds  Number".  Their 
model,  suitably  modified,  was  used  for  the  present  analysis. 

The  response  of  the  base  bending  moment  of  a  slender,  vertically- 
mounted  cantilever  cylinder  to  a  turbulent  flow  field  with  shear  was  measured 
for  the  drag  (streamwise)  and  lift  directions,  and  the  drag  response  was  com¬ 
pared  to  the  response  calculated  for  .  same  flow  field  ty  a  linear  theory.-due 
to  B .  Etkin. 

This  theory  makes  use  of  the  method  of  normal  modes,  and  includes  only 
the  effect  of  random  turbulence  in  the  generalized  force  function.  The  turbu¬ 
lent  field  was  determined  by  measurement  of  the  requisite  space-time  corre¬ 
lations  . 


Agreement  between  the  theoretical  predictions  and  experimental  meas¬ 
urements  was  found  to  be  quite  good,  especially  at  lower  frequencies. 
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INTRODUCTION: 


Considerable  interest  has  been  shown;,  Tr.  rece.n*.  yea’-s.  ir.  the  v espor.s 
of  tall,  slenaer  structures  (radio  antennae  ,  tall  towers,  smoke  stacks,  launch 
vehicles,  etc)  to  random  excitation  by  atmospheric  turbulence 

The  response  to  any  such  excitation  is  also  random  and  so  must  be 
treated  in  a  statistical,  manner.  This  response  can  be  described  in  terms  of 
either  a  spectral  density  functicn  (in  the  frequency  domain)  or  a  correlation 
function  (in  the  time  domain).  Following  usual  practice,  she  response  analysis 
is  expressed  in  the  frequency  domain. 

This  report  compares  the  computed  spectral  density  with  the  measured 
spectral  density  of  the  drag  response  for  a  thin-walled  uniform  brass  cylinder 
which  has  mechanical  and  aerodynamic  properties  similar  to  the  contemporary 
slender  structures  mentioned  above . 

The  turbulent  stream  can  be  looked  upon,  as  consisting  of  two  compo¬ 
nents,  one  static  and  one  dynamic.  That  is,  we  have  an  os  :  cl hating  turbulent 
velocity  component  superimposed  on  some  mean  flow,  in  this  case,  to  provide 
a  better  test  of  the  mathematical  model  proposed,  a  nor.- uniform  mean  velocity 
profile  (shear  profile)  was  employed.  The  shear  turbulent  flow  field  was 
developed  behind  a  "shear"  grid  (see  Fig.l)  and  its  mean  velocity  profile  is 
shown  in  Fig. 2,  while  the  associated  turbulent  component  is  recorded  in.  Fig, 3- 

The  turbulence  properties  were  initially  determined  in  the  time  domain, 
because  the  availability  of  a  correlation  function  computer  (FAR  MODEL  100) 
made  these  measurements  most  expedient.  A  digital  Fourier  transformation  of 
these  correlations  was  employed  to  generate  the  c.o'  responding  spectrum  function 

Just  as  the  flow  field  consists  of  a  static  and  dynamic  component  so 
the  bending  response  of  the  beam  in  the  d~«.g  direction  has  a  mean  and  a  time 
varying  component.  It  is  this  fluctuating  drag  component  that  is  of  primary 
interest. 

1.  EXPERIMENT  AND  ANAL  Yd  .13 
1 ■ 1  Model 

The  thin-walled  brass  cylinder  (sec  Fig.4)used  for  the  model  originally 
had  the  following  dimensions  (ex  Campbell  model,  reference  3) « 

Length  29 „ 5  inches 

Diameter  ‘•^he. 

Wall  thickness  010  inches 

Initial  measurements,  however,  showed  that  the  structural  damping  was 
lower  than  wanted  (£„  -  .001),  so  an  outer:  sheath  of  visco-elastic  material 

was  heat  shrunk  onto'"5*  the  cylinder  (see  Fig  4).  Th.;s  modification  more  than 
doubled  the  structural  damping  ( L  “  .0022),  yet  did  not  increase  bending 
stiffness  measurably.  The  final  dimensions  cf  the  mode,  were 


Length 


29. r:  inches 


Diameter 
Wall  thickness 


.‘>0  inches 
.026  inches 


1.2  Instrumentation 


The  principal  instrumentation  used  to  obtain  the  bending  response  data 
was  strain  gauges.  Two  four-gauge  bridges,  with  one  gauge  per  bridge  arm, 
were  mounted  on  the  model  near  the  base  (see  Fig. 5).  The  gauges  were  powered 
by  a  12  volt  d.c.  source  and  their  outputs  were  multiplied  by  a  gain  of  500 
in  highly  accurate  differential  amplifiers.  The  output  from  20  cps  to  500 
cps,  the  upper  portion  of  the  range  of  interest,  was  analysed  on  line.  The 
lift  signal  (with  zero  mean)  was  frequency  analyzed  directly  by  a  Bruel  and 
Kjaer  wave  analyzer  and  random  noise  meter,  while  the  drag  signal  “irst  had 
the  mean  signal  removed  (and  recorded)  before  frequency  analysis. 

The  wave  analyzer  was  of  the  constant  percent  bandwidth  type,  and  the 
bandwidth  used  was  the  most  narrow  (7.5  percent).  The  random  noise  meter  was 
equipped  with  a  variable  time-constant  and  was  used  to  obtain  the  true  rms 
signal  from  the  wave  analyzer,  whose  low  frequency  capability  extended  to 
20  cps. 


Because  analysis  was  required  down  to  as  low  a  frequency  as  possible, 
some  data  was  tape  recorded  at  1-7/8  ips  on  an  Ampex  SP  300  recorder .  By 
playing  this  record  back  at  15  ips,  analysis  was  possible  down  to  2-5  cps. 

In  this  manner  a  frequency  range  front  2.5  to  500  cps  was  attained. 

1.3  Mount  and  Installation 

Care  was  required  in  the  manner  of  supporting  the  model  so  that  ex¬ 
traneous  vibrations  in  the  frequency  band  of  interest  would  be  eliminated  . 
Campbell' 3  mount  was  used  as  it  was  entirely  adequate. 

A  30  pound  weight  was  supported  by  four  thin  steel  rods  (l/8  inch 
diameter)  as  shown  in  Fig. 6,  and  the  cylinder  was  fastened  to  the  weight. 

The  natural  frequency  of  the  mount  was  i.3  cps. 

The  mount  was  then  isolated  from  the  wind  tunnel  by  being  supported 
directly  from  the  concrete  floor  beneath  the  tunnel  (see  Fig’s.  6  and  7), 

A  large  fairing  with  an  attached  ground  board  was  used  to  shield  the 
weight  from  the  flow,  and  to  remove  the  effects  of  the  tunnel  boundary  layer. 
The  fairing  was  10  arranged  that  the  lower  portion  of  the  cylinder  on  which 
the  strain  gauges  were  mounted  was  also  shielded,  so  that  exactly  27*0  inches 
of  the  cylinder  extended  above  the  ground  board  (see  Fig. 8). 

The  cylinder  was  located  in  the  diffuser  section  downstream  of  the 
test  section,  125  inches  from  the  shear  grid.  The  tunnel  aerodynamic  outline 
is  shown  in  Fig. 9  &nd  the  grid  and  test  location  are  both  indicated. 

1.4  Scaling 

The  requirements  for  dynamic  and  geometric  similarity  of  an  aeroelastic 
model,  from  dimensional  analysis,  were  found  to  be 


total  damping  ratio 


Et/p  U2 


m/p  bc 


elasticity  parameter 
density  parameter 


pU  b/p  Reynold's  number 


Flow  structure 


b jl  Geometric  similarity 

where  the  model  and  full  size  structure  must  have  the  same  values  of  the 
above  parameters. 

Assuming  the  same  mean  velocity  profile  for  both  model  and  full  scale 
(no  attempt  was  made  to  scale  the  earth's  boundary  layer)  then  the  model 
represents  one  possible  steel  structure  of  the  following  dimensions  - 


length 

diameter 

wall  thickness 


216  feet 
4.0  feet 
•  5  inches 


1.5  Theoretical  Bending  Moment  Spectrum  (after  Campbell  &  Etkin) 


Appendix  B  gives  the  general  theory,  and  by  using  some  simplifying 
assumptions  the  theory  was  applied  to  the  cylindrical  model. 

The  basic  deflection  spectral  relation  is  given  by  (b2l)  as 
00  00 

<t>  (z,if)=£  Z  F  (z)F  (z)G*  (if)  G(if)<fr_  „  (f)  (1-5-1) 

xx  m=l  n=l  m  n  m  n  3m3m 

where  Fm(z),  F^z)  are  the  raotie  shaPe  functions;  Gm(if),  Gn(if)  are  the 
transfer  functions  (*  indicates  the  complex  conjugate)  for  the  rath  end  nth 
modes  and  $3^^?)  ts  the  generalized  input  force  cros.  -spectrum  function. 

For  the  case  of  a  lightly  damped  system  (which  we  have),  and  for  a 
structure  whose  natural  mode  frequencies  are  well  separated  we  can  assume 
that  there  is  no  cross  coupling  between  modes  of  the  model.  Thus 

%  ,  (?)  =  0 

^  m^m 

for  m  f  n.  Hence  the  above  equation  reduces  to 


<*>  (z,if)  =  £  F  Hz)  G  (if)  2  <*> 


■^n^n 


(1-5.2) 


The  bending  moment-deflection  relation  is 
M  (z,*)  =  B(z)  d2  x/o  z2 


(1.5.3) 


3 


where  B(z)  =  El(z)  =  constant.  But  from  the  modal  representation  of  deflection 


Hence 


x(z,t)  -  £  F  (z)  q  (t) 


£2x/dz2  =  £ 


n  ‘  '  n 


.  2 


n  dz 


qn(t) 


(1.5,4) 


i.e.,  cJ2x/dz2  =  £  F  "  (z)  qn  (t) 

n  n 


For  the  case  of  base  bending  moment 


ff"  ■  1  V  <•)  «„<*> 

oz  z  =  0  n 


(1.5*5) 


expressed 


Using  (1.5.3)  and  (1.5.5)  the  bending  moment  spectrum  may  be 


Vf>  ‘  §  ^  Pn"2<°>  |  °n<lf>  T  *3J  <f> 

JnJn 


(1.5.6) 


The  values  of  B,  Fn"  (o)  and  |  G^if)  |  are  given  in  Appendix  C. 

The  generalized  force  spectrum  <£-  g  (f)  is  given  by  equation  (b  25)  and  is 
for  this  case  n  n 

V 

<t>3  (f)  =  JJ  0uu(a,p,f)  jgn(a)gn(P)  +  4tt2  f\(<*)hnO) 

o 

+  i  2  Trf  [gn(a)hn(p)-gnO)hn(a)]  j*  da  dp  (1.5*7) 

Now  consider  the  integrals  over  the  (a,p;  domain  illustrated  in 
Fig. 10.  Let  the  integrals  be  evaluated  by  summing  at  pairs  of  elements  (p,q) 
that  are  symmetric  with  respect  to  the  diagonal.  Since  (a, (3)  are  interchanged 
at  p  and  q,  the  following  relations  hold 


h  (a)  h  (p) 
n  '  n'  ] 


en(=)  g„(f' 


hn(a)  hnO) 


(1.5.8) 


Using  these  relations  plus  the  symmetry  relations  (a  l6)  and  (a  17) 
and  equation  (a  11)  (which  separates  0uu(a,P,f)  into  real  and  imaginary  parts) 
in  equation  (l.5.7)^ve  obtain  the  generalized  force  spectrum  function 

<t>3  ^  (f)  =  j  j  |  01(a,p,f)[gn(a)gn(p)+  4/f2hn(a)hn(p)  ] 

+  27rf  02(a,p,f)[gn(a)hn(p)-gn(p)hn(a)jj  da  dp  (1.5. 9) 


4 


i 


Where  l'  is  the  exposed  length  \ 27.0  inches)  of  the  cylinder  and  0^  and  0g 
are  the  real  and  imaginary  components  of  the  velocity  spectrum  function.  An 
estimate  of  1 h ;  true  velocity  cross  spectrum  is  required  to  evaluate  (1.5.9). 

At  this  point  two  assumptions  which  have  been  important  in  the 
development  of  the  theory  should  be  emphasized. 

(i)  The  fetructure  is  sufficiently  slender  for  the  secondary  span- 
wise  flow  and  consequent  redistribution  of  pressures  to  be 
neglected,  such  that  the  pressures  on  any  section  of  the  span 
are  only  due  to  the  wind  incident  on  that  section.  This  is 
the  "strip  theory"  assumption. 

(ii)  The  turbulent  fluctuations  are  so  small  compared  to  the  mean 
velocity  that  the  fluctuating  loads  can  be  expressed  as 
linear  functions  of  the  turbulent  velocities. 


Both  of  these  will  be  discussed  in  section  2.  Assuming  two-dimen¬ 
sional  flow  over  a  large  portion  of  such  a  long  slender  cylinder  the  drag 
coefficient  chosen  was  <'Ref.l5) 

CD  =  1*2 

1.6  Theoretical  Mean  Bending  Moment 


The  force  per  unit  length  on  the  cylinder  due  to  the  mean  flow 
at  position  z  is 

-  Cl>Pb  U2  (z)  •'  c 


f  (z)  = 


(1.6.1) 


Thus  the  force  on  the  increment  dz  at  z  is 

C  pb 

dF  =  f(z)  dz  =  — 5 - 


U  2(z)  dz 


(1.6.2) 


The  bending  moment  due  to  the  force  on  the  increment  dz  at  z  is 
_  C  Pb 

dM  =  dF  Z  =  IT  (z)  Z  dz 


(1*6.3) 


Hence  the  total  mean  bending  moment,  is 

=  _  vb  r29'5*, 


f  f  «■<•> 


(1.6.10 


The  limits  are  due  to  the  fact  that  not  all  of  the  cylinder  is  exposed,  so 
U  below  z  =  2.5  inches  is  zero. 
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1.7  Experimental  Bending  Moment  Spectra 


The  experimental  bending  moment  spectra  can  be  found  from  the  rms 
output  of  the  wave  analyser.  The  filter  characteristic  of  the  wave  analyzer 
for  sine  wave  inputs  at  frequency  f  is  shown  in  Fig. 11. 

mean  square  output  e  2 

g  (f/fj  =  -  -  (1-7-1) 

mean  square  input  2 

ei 

where  fs  ia  the  centre  frequency.  Then  the  mean  square  output  signal  from  the 
wave  analyzer  is 

e02  (f,  fs)  =  g(f/fs)  ei2(f)  (1-7.2) 

Now  the  strain  gauge  aignalj  e^  (see  Fig. 67),  is  directly  proportional  to  the 
base  bending  moment,  i.e., 

ei  =  ri  M  (1.7-3) 

and 

et2  =  n2  M2  (1.7 .4) 

where  tj  is  the  slope  of  the  strain  gauge  calibration  curves  (see  Fig. 12).  Thus 
the  wave  analyzer  output  signal  is 

eQ2  (f,fs)  =  g(f/fsh  2  M2  (1.7.5) 

But  from  the  bending  moment  spectrum  curve  the  mean  square  bending  moment  about 
the  mean,  for  a  wide-band  input  is 

00 

^  =  2  J  (f)  df  (1.7.6) 

Thua  the  output  corresponding  to  the  wave-analyzer  centre  frequency  f  is 

s 

00 

eQ2  (fg)  =  2  7]2  f  0BM(f)  g(f/fs)  df  (1.7.7) 

Assuming  0jjM(f‘)  is  approximately  constant  over  the  frequency  band  of  the  filter, 
(1.7.7)  becomes 


00 


«02  <fs>  -  2  'i2  0BM(fs)/  aW*B)  « 

w  0 

(1.7.8) 

00 

Jg(f/fs)  d(f/fa)  =  G 

(1.7.9) 

O 
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Then 


2  Vf) = 


(1.7.10) 


e  2  (f  )  /  rf 
o  's'  1 


f  G 
s 


For  the  filter  used  G  -  .0763. 

The  mean  bending  moment  is  obtained  from  the  mean  voltage  output 
as  developed  in  the  analogue  analysis. 

2.  TURBULENCE  DATA 

2.1  Measurement  of  Ruu 


The  turbulent  flow  field  was  developed  behind  a  shear  grid  (Fig.l) 
in  the  OTIAS  Subsonic  Wind  Tunnel. 

The  initial  survey  was  carried  out  using  a  single  hot  wire  probe, 
ranged  in  increments  of  1  inch  from  the  ground  board  (z  =  .125  in)  to  z  =  30  in. 
Employing  the  analogue  computation  scheme  shown  in  Fig. 13,  (mechanized  on  the 
UTIAS  PACE  computer)  the  mean  was  extracted  from  the  signal  and  displayed  on  a 
digital  voltmeter,  while  the  rms  value  of  the  fluctuating  component  was  measured 
by  a  Bruel  and  KJaer  random  noise  meter.  The  linearized  hot  wire  output  was 
converted  to  velocity  using  the  calibration  curve,  Fig.l^.  The  resulting  mean 
velocity  and  turbulence  intensity  profiles  are  shown  in  Fig’s.  2  and  3  res¬ 
pectively.  (Figure  27  shows  Reynolds  number). 

The  turbulence  intensity  over  the  first  6-8  inches  of  the  cylinder 
seems  too  high  for  the  linear  assumption  (ii)  of  (l„5)  to  hold.  However,  a 
constant  temperature  anemometer  tends  to  over  estimate  the  intensity  (see 
Ref .t)  of  high  intensity  turbulence  so  the  real  case  is  not  as  bad  as  first 
appears.  Also  this  phenomenon  occurs  over  the  lower  portion  of  the  cylinder 
where  the  drag  and  contribution  to  the  bending  moment  are  small  so  the  error 
Incurred  would  be  small  if  assumption  (ii;  is  retained. 

One  other  result  of  such  high  intensities  is  that  reverse  flow, 
which  the  probe  cannot  separate  from  forward  flow,  is  possible  part  of  the 
time,  but  this  effect  is  small  also.  For  example  (assuming  a  normal  distri¬ 
bution  for  u(x,t)  at  z  =  3  inches,  the  probability  of  being  within  2cr  =  2u/U=  1.0 
is  95 .*+$.  Thus  reverse  flow  occurs  cnly  2.3$  of  the  time.  At  z  =  8  inches  the 
probability  of  reverse  flow  is  oniy  .15$  and  at  z  =  .125  inches  the  probability 
is  21.5$.  Obviously  this  has  little  or  no  effect  on  the  response. 

For  the  cross  correlation  measurements  (see  Appendix  Al),  two 
probes  were  required  (see  Fig. 15).  The  signals  from  these  probes  first  had  the 
mean  components  removed  (as  above)  and  were  then  fed  into  a  Princeton  Applied 
Research  Correlation  Function  Computer  (Model  100) .  The  probe  arrangement  and 
notation  used, is  diagrammed  below. 


r 


Probe  Q  is  the  fixed  probe,  while  probe  (3  ranges  around  each  Q  posi¬ 
tion.  Increments  in  Q  and  p  of  2.7  inches  were  used  for  a  survey  range  of 
.125  (as  close  to  the  ground  board  as  possible)  to  27.0  inches  (i.e.,  .125 
2.7,  5-4 . 27.0). 


The  PAR  computer  calculates  the  correlation  function  for  100  points 
at  lag  increments  of  TmHY/lOO,  where  Tnjax  is  the  maximum  lag  time  chosen.  This 
choice  is  based  on  both  required  frequency  response,  and  the  desire  to  include 
the  whole  correlation  curve.  A  value  of  =  .2  sec  (At  =  .002  sec)  was 

chosen  as  this  provided  a  frequency  cutoff  at  250  cps.  which  covered  the  first 
three  bending  modes,  and  included  the  majority  of  the  correlation  curve.  The 
averaging  time  for  the  computer  was  4o  seconds.  The  correlation  function 
computed  may  be  written 


fuu  (a.?,  t) 


U 


40 


u  (a,t)  u  (« 


►-» 


n)  dt 


^  -1  1  \ 

r-  .  1 .  >.  / 


where  n  =  1,2, - 100  (i.e,  t  =  0,  .002, - - .2  sec).  As  can  be  seen 

from  equation  (a  18)  this  is  only  one  half  of  the  total  cross  correlation.  By 
reversing  the  probe  inputs  to  the  PAR  computer  we  get  from  (2.1.1). 

l  r 40 

Ruu(P»  Q’T>  =  fe  /  u  u<a’  t  +  .002  n  )  dt  (2.1.2) 

J  o 


However,  applying  the  symmetry  properties  (a  i.4)  and  (a  15)  to 
(2.1.2)  we  see  that 
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(5.1.3) 


(p,  a,  t)  =  R  (q,P,-t) 

Thus  by  reversing  the  inputs,  the  second,  half  of  the  correlation  curve  was  ob¬ 
tained. 

These  symmetry  properties  were  also  used  to  reduce  the  number  of 
measurements  required.  Fran  (2.1. 3)  we  can  see  that  reversing  the  time  scale 
for  a  fixed  (a,  fi)  relationship  gives  the  cross  correlation  for  the  (P,o0  in¬ 
put.  Instead  of  having  to  measure  the  entire  tensor  shown  below,  we  only  need 
measure  the  shaded  half. 

- - - •*- 


Figures  16  to  26  show  a  typical  set  of  correlation  curves,  obtained 
for  a  =  13.5  inches  and  P  ranging  from  .125  inches  to  27. 0  inches  (the  .125 
inch  value  is  considered  zero) . 

2.2  Turbulence  Spectra 

The  velocity  cross  correlations  R*uu(0!,P,t)  were  Fourier  transformed 
as  in  equations  (a  12)  and  (a  13)  and  modified  (Hanning  factor,  equations  (a  2b) 
applied)  digitally  as  described  in  Appendix  A. 2.  The  transformed  spectra  corr¬ 
esponding  to  the  correlations  shown  in  figures  16  to  26  are  shown  in  figures 

28  to  38. 

The  total  array  of  velocity  cross  spectra  was  employed  with  equation 
(1.5*9)  to  calculate  the  generalized  input  force  spectra  for  each  particular 
mode.  Simpsoris  rule  integration  was  employed  in  the  digital  program  and  should 
provide  sufficient  accuracy.  The  resulting  inputs  can  be  seen  in  figure  39* 

They  demonstrated  the  form  to  be  expected.  The  first  mode  input  has  the  greatest 
content  at  low  frequencies  and  dropa  off  most  rapidly  with  frequency.  The 
second  mode  input  is  intermediate  in  magnitude  and  frequency  range,  while  the 
third  mode  forcing  function  has  the  lowest  content  at  low  frequencies,  but  has 
the  largest  high  frequency  components. 

Unfortunately  all  the  signals  disappeared  at  the  noise  level  around 
90  to  125  cps,  so  the  third  mode  peak  for  the  bending  moment  spectrum  could  not 
be  calculated. 

Also  due  to  the  finite  maximum  lag  time  available,  some  of  each 
correlation  curve  was  lost,  and  in  spite  of  the  Hannirg  correction,  the  low 
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frequency  values  will  consequently  be  somewhat  low. 

Some  workers  (references  ?  and  6)  have  used  a  coherence  function 
(normalized  co-spectrum)  to  calculate  the  cross  spectral  densities  given  the 
spatial  distribution  of  velocity  power  spectral  densities.  They  have  assumed 
that 

R  =  0L(a,p,f)/  [0uu(a,a,f)  0uuO.P»f)3  *  (2.2.1) 


when  plotted  versus 

Up)  (2.2.2) 


depends  on  the  turbulent  field. 

Vickery  (reference  6)  quotes  a  value  of  k  =  9  for  rolling  wooded 
country,  and  Davenport  (reference  5)  quotes  k  =  6  for  wooded  countryside 
(after  Singer)  and  7.7  for  open  grassland. 

Figure  ^0  presents  some  of  the  turbulence  data  for  various  probe 
stations  plotted  in  the  form  given  by  equations  (2.2,1)  and  (2.2.2)  .  The 
curve  fitted  through  the  points  closely  approximates 

E  =  e’6-657 


f  =  r<wi(V 

may  be  expressed  in  the  form 


where  k  is  some  constant  value  which 


This  approximation  would  seem  of  considerable  value  as  it  greatly 
reduces  the  work  required  to  estimate  the  cross  spectra. 


The  power  spectra  could  be  measured  directly  and  it  is  felt  that 
any  errors  incurred  in  the  exponential  approximation  (2.2.3)  would  be  outweighed 
by  the  very  good  accuracy  achievable  in  the  measurement  of  the  power  spectra, 
as  compared  to  first  obtaining  the  correlation,  and  then  transforming  to  the 
frequency  domain. 


2.3  Scale  of  Turbulence  (see  Appendix  A. 3) 

The  ratios *R  (a,p,0) 

zr- -  (2.3.1) 

RUU 

which  represent  the  normalized  cross  correlations  are  plotted  versus  probe 
separation  (a-(3)  for  all  eleven  survey  positions  in  Figs,  bj.  to  51.  The  areas 
under  these  curves.’ ;epresent  the  lateral  integral  scales  L  given  by  equation 
(a  28).  The  distribution  of  scale  over  the  cylinder  length  is  shown  in  Fig.  52. 


Taylor’s 


The  longitudinal  scale  given  by  equation  (a  28)  is  computed  (assuming 
hypothesis)  by  finding  the  area  under  the  autocorrelation  curve 

A(a)  =  j 
u  o 


Ku  T) 


dT 


(2.3.2) 
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and  then  multiplying  by 


u  (a)  /^uu  (a,a,o) 

to  normalize  the  correlation  and  convert  to  length  units  so  that  the  result. 

A  (a)  U  (a) 

R*  (a,a,o) 

corresponds  to  Lx(z)  (equation  (a  27)).  These  are  plotted  versus  z  in  Fig. 53. 

For  isotropic  flow,  Lx(z)  =  2  Lz(z)  (reference  10).  This  does  not 
hold  for  the  present  anisotropic  flow,  as  can  be  seen  in  Fig's.  52  and  53. 
Generally  however,  Lx(z  )>  Lz(z)  except  at  the  base  of  the  cylinder  where  the 
turbulence  intensity  is  very  high,  and  here  L„(z)  >  Lv(z). 

For  similarity  between  model  and  full  scale  flows,  the  ratio 

r  Lx  (z)  i 

should  be  the  same  for  model  and  structure. 

Since  no  attempt  was  made  to  scale  flows,  this  parameter  has  no 
application  here  other  than  to  support  the  strip  theory  assumption  which 
requires  L  »  b. 

3  •  EXPERIMENTAL  RESULTS 


3.1  Mean  Bending  Moment 

The  mean  bending  moment  for  the  lift  direction  was  zero,  >.nd  that  of 
the  drag  direction  was  found  to  be 

Mp  -  3.20  in  lb. 


3.2  Spectral  Response 

The  response  spectra  were  calculated  as  outlined  in  section  1.7.  The 
lift  response  is  shown  in  Fig, 55.  The  magnitude  of  the  third  mode  peak  is 
significantly  greater  than  that  of  the  second  mode  due  to  vortex  shedding.  The 
variation  with  height  of  Strouhal  shedding  frequency,  assuming  that  a  Strouhal 
number  of  .2  is  applicable  to  the  situation,  is  shown  in  Fig.51*.  As  can  be 
seen,  over  the  top  half  of  the  Cylinder  the  shedding  frequency  ranges  from  150 
cps.  to  250  cps.  The  third  mode  natural  frequency  is  213  cps.,  so  this  mode 
is  very  likely  being  excited  by  the  vortex  shedding.  The  variance  of  the  lift 
signal  was  equivalent  to 

o  2  =  1.18  (in  lb)2 

L 

The  drag  response  spectrum  is  given  in  Fig. 56.  As  expected,  the  drag 
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response  is  greater  than  that  in  the  lift  direction,  in  spite  of  the  vortex 
excitation.  The  variance  for  the  drag  signal  was  equivalent  to 


u2  =  6.17  (in  lb)2 


The  ratios  of  peak  values  of  the  drag  spectrum  for  the  va.  ious  inodes 
are  (order  of  magnitude  only) 


first  to  second 

10' 

first  to  third 

J 

10 

second  to  third 

10 

For  the  cases  where  vortex  excitation  was  present  Campbell  obtained 

first  to  second  10^ 

first  to  third  10^ 

second  to  third  10 ^ 

When  no  vertex  excitation  was  present  the  ratio  "second  to  third"  was  approxi¬ 
mately  10  . 

This  would  suggest  that  as  in  the  previous  report,  there  is  some  cross¬ 
coupling  between  lift  and  drag  in  the  third  mode  response. 

4.  THEORETICAL  RESULTS 


4.1  Mean  Bending  Moment 

The  mean  bending  moment  was  computed  as  shown  in  section  1.6.,  using 
a  graphical  evaluation  of  equation  (l.6.4).  The  value  obtained  was 

Mp  =  3-38  in  lb. 

which  is  5.8  percent  higher  than  the  measured  value.  The  bending  moment  calcula¬ 
tion  was  performed  for  a  structure  perpendicular  to  the  mean  flow.  Actually 
there  would  be  an  equilibrium  between  the  aerodynamic  load  and  elastic  restoring 
force,  in  which  the  structure  would  be  deflected,  and  so  the  load  on  the  curved 
cylinder  would  be  slightly  reduced  due  to  the  effective  lowering  of  the  normal 
flow  velocity  over  the  cylinder.  (The  magnitude  of  this  effect  was  not,  however, 
calculated) . 

4.2  Spectral  Response 

A  digital  program  was  written  to  evaluate  equation  (1.5.6)  for  (f) 
using  the  generalized  input  force  spectra  calculated  from  equation  (1.5. 9). 

The  computed  drag  response  is  shown  in  Fig. 57-  However,  this  spectrum 
as  it  stands,  was  not  felt  to  be  sufficiently  useful. 

An  assumption  was  made  in  the  calculation  of  the  experimental  spectrum, 
that  0gM(f)  was  constant  over  the  bandwidth  of  the  wave  analyzer.  Obviously, 
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this  could  not  be  true  near  a  narrow  response  peak  for  a  7-5  percent  bandwidth. 
Since  the  wave  analyzer  has  a  finite  bandwidth,  the  spectrum  at  a  particular 
frequency  point  is  really  an  average  over  a  frequency  band  around  that  frequency. 
This  point  would  tend  to  confuse  any  direct  comparison  between  theory  and  experi¬ 
ment  for  it  would  reduce  experimental  peak  values,  and  increase  the  values  between 
the  peaks.  As  the  point  of  this  paper  is  just  such  a  comparison,  it  was  felt 
necessary  to  filter  the  computed  data  digitally,  using  the  filter  shape  of  the 
analyzer  (Fig. 11).  Then  the  data  would  be  directly  comparable,  and  any  discre¬ 
pancies  would  not  be  due  to  the  measurement  and  analysis  technique  of  the 
experimentally  obtained  spectrum.  The  filtered  theoretical  drag  spectrum  is 
compared  to  the  unfiltered  spectrum  in  Fig. 57* 

One  probable  cause  of  the  loss  of  the  force  spectrum  above  100  cps. 
was  an  insufficiently  high  signal  to  noise  ratio  in  the  measurement  of  the 
turbulence.  At  the  same  time,  other  workers  at  the  Institute  have  found  that 
the  validity  of  spectra  obtained  from  the  correlation  function  by  Fourier 
transformation  were  unreliable  once  the  spectrum  v*>s  dcwn  from  the  peak  value 
by  three  orders  of  magnitude.  This  is  one  of  tin  reasons  for  the  suggestion 
(section  2.2),  that  the  exponential  relation  between  power  spectra  and  co¬ 
spectra  be  employed. 

The  variance  for  the  frequency  range  of  2  cps  to  100  cps.  is 

(  %  )  =  6.43  (in  lb)2 

F 

where  the  T  superscript  represents  a  "truncated"  spectrum  from  2  to  100  cps, 
and  the  subscript  F  represents  the  filtered  computed  response  (see  next 
section) . 

5.  COMPARISON  OF  THEORY  AMD  EXPERIMENT 


Figure  56  directly  compares  the  filtered  theoretical  and  experimental 
drag  bending  moment  spectra.  The  comparison  is  facilitated  by  reference  to 
Fig. 58,  whi^h  gives  a  plot  of  the  ratio  of  the  filtered  theoretical  and  experi¬ 
mental  spectra.  We  can  see  that  the  agreement  is  reasonable  (i.e.,  ratio 
reasonably  close  to  1.0)  up  to  about  20  cps.  but  is  poor  at  the  highest  fre¬ 
quencies.  Surprisingly  the  theoretical  values  are  lower  than  the  experimental 
ones  (one  might  expect  larger  values)  except  in  the  frequency  ranges  near  the 
peaks,  wh«  e  the  theoretical  values  are  too  high.  This  immediatly  suggests 
that  the  estimates  of  the  generalized  forces  and  the  damping  might  have  been 
in  error. 


Underestimation  of  the  generalized  force  spectrum  would  lower  the 
bending  moment  spectrum  over  the  entire  frequency  range  and  partially  account 
for  the  low  portions  between  the  peaks  of  Fig. 56.  However,  near  a  peak,  the 
damping  term  in  the  transfer  function 

’  1/6  hv  f2  f  2  (  £a  2  +  5s  2) 

‘  n  '  s  n  5  n  ' 


dominates  as  f  approaches  f^. 

If  the  damping  estimate  were  low,  it  would  cause  little  change  away 
from  the  peaks,  but  at  a  peak  could  significantly  raise  the  peak  value 
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The  structural  damping  was  felt  to  be  quite  weil  delict"!,  so  any  error 
in  the  damping  is  most  likely  due  to  the  estimation  of  the  aerodynamic  part. 

Both  the  input  forte  spectrum  and  the  aerodynamic  damping  were  bused  on  the 
initial  assumptions  of  strip  theory  and  linearization.  Certainly  the  lineari¬ 
zation  fails  to  some  extent  as  does  the  two  dimensional,  assumption,  and  so  one 
would  expect  r.i  erroneous  theoretical  estimate,  though  to  the  nigh  side. 

The  agreement  appears  to  fail  progressively  with  increasing  frequency. 
It  is  possible  that  the  interaction  of  the  flow  with  the  more  complex  second 
mode  shape  causes  some  of  the  discrepancy.  This  could  be  readily  determined  by 
utilizing  a  stiffer  cylinder  whose  first  mode  natural  frequency  would  be  in 
the  neighbourhood  of  75  cps. 

It  has  been  suggested  above  that  - 

(i)  the  generalized  force  spectrum  computed  is  too  small  in  magnitude, 
especially  for  the  second  mode  input. 

(ii)  the  aerodynamic  damping  is  too  small  also. 

Both  these  quantities  are  dependent  on  one  very  important  parameter, 
the  drag  coefficient  Cq.  The  value  used  was  1.2,  which  is  typical  of  the  flow 
Reynolds  numbers  (see  Fig.?7)  for  the  uniform  two  dimensional  flow  expected  over 
a  large  portion  of  such  a  long  slender  cylinder . 

In  a  turbulent  field  a  somewhat  lower  value  might  be  expected.  However 
Scruton  (Ref. 9)  has  shown  that  the  drag  coefficient  of  a  !or:g  slender  finite 
cylinder  increases  quite  sharply,  at  least  25-30  percent,  toward  the  tip.  On 
the  basis  of  this  information,  it  is  very  possible  that  the  drag  coefficient 
has  been  underestimated,  and  as  this  has  very  important  consequences,  should  be 
investigated  more  closely. 

An  increase  in  would,  of  course,  increase  the  generalized  input 
force  spectrum  and  so  improve  the  agreement  between  theory  and  experiment  away 
from  the  peaks.  At  the  peaks,  an  increase  in  Cp  would  cause  a  decrease  in  the 
theoretically  computed  values.  This  occurs  because,  although  the  input  forces 
are  increased,  the  larger  leads  to  a  relatively  greater  reduction  in  the 
transfer  function,  through  the  aerodynamic  damping  terms  (see  equation  (5.1)). 

For  example,  the  ratios  of  the  peak  values  for  0  -  Lo  to  the  peak 

values  for  =  1.2  for  the  first  3  modes  are 

Mode  Ratio 

1  .875 

2  .950 

3  .960 

As  might  be  expected,  the  Cp  change  has  a  smaller  effect  for  the  higher 
modes  because  the  aerodynamic  damping  is  a  smaller  and  smaller  component  of 
the  total  damping. 

Good  agreement  up  to  50  cps.  would  require  a  value  of  C-,  of  about.  1.45 
to  1.50.  However,  no  reasonable  adjustment  in  Cq  would  correct  )Sr  the  total 
error  at  the  second  mode  peak. 
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The  data  presented  beyond  90  cps.  was  expected  to  be  low,  as  the  in¬ 
put  force  spectra  were  approaching  noise  level. 

Table  3  provides  a  comparison  of  the  contribution  to  the  total  variance, 
of  various  frequency  bands  of  the  filtered  and  unfiltered  theoretical  response, 
with  the  experimental  response.  The  effect  of  filtering  is  seen  as  a  small 
increase  in  the  variance  components  resulting  in  a  slightly  higher  filtered 
variance.  (4$).  These  ratios  were  plotted  in  Fig.  59  and,  of  course,  reflect 
the  form  of  Fig. 56. 


The  total  variance  for 
was  found  to  be 

T 

[  "n  J  = 


the  theoretical  truncated  filtered  spectrum 
6.43  (in  lb)2 


and  the  variance  for  the  truncated  experimental  spectrum  was  found  to  be 

T 

[  <£  3  =  5-99  (in  lb)2 

EXP 

Vickery,  reference  6,  found  that  his  computation  for  a  somewhat 
similar  case  resulted  in  the  theoretical  variance  being  high  by  about  22.%. 

In  this  case  the  theory  is  %  high. 

Again  it  should  be  pointed  out  that  the  discrepancy  in  the  region  of 
the  second  mode  has  not  been  explained,  and  it  is  suggested  that  this  point  be 
the  subject  of  further  research. 

6.0  CONCLUSIONS 

(i)  The  theory  provides  a  good  description  of  the  response  of  a 
slender  structure  to  turbulence  at  the  first  mode  for  suffi¬ 
ciently  low  frequencies  of  vibration.  The  variance  of  the 
theoretically  calculated  response  was  eight  percent  higher 
than  the  measured  value. 

(ii)  At  higher  modes  and/or  for  sufficiently  high  frequencies  the 

theoretically  predicted  values  of  drag  response  are  low.  This 
discrepancy  may  be  caused  by  some  or  all  of  - 

a)  Failure  of  the  two-dimensional  assumptions  at  higher  frequencies 

b)  A  low  value  of  C^,  although  no  reasonable  value  could  explain 
the  total  discrepancy. 

c)  Some  flow-response  phenomenon  of  the  higher  modes  not  represented 
by  the  theory. 

(iii)  The  two-dimensional  subcritical  value  assumed  for  the  drag  co¬ 
efficient,  Cj  =  1.2,  appears  to  be  low. 

(iv)  The  overestimation  of  the  response  at  the  natural  mode  peaks  may 
be  due  to  a  low  estimate  of  aerodynamic  damping  which  can  be  tied 
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to  a  low  value  of  C^.  However,  this  does  not  explain  the 
entire  response  at  the  second  mode  peak,  and  so  it  is 
possible  that  the  value  of  structural  damping  for  the  second 
mode  is  in  error. 

(v)  The  lift  response  to  turbulence  was  shown  to  be  considerably 
lower  than  the  drag  response. 

(vi)  More  work  is  required  to  determine  the  proper  value  (s)  for 

Cjj  and  to  investigate  the  high  frequency,  higher  mode  response. 
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APPENDIX  A: 


STATISTICAL  SPECIFICATION  OF  THE  TURBULENT  FIELD 


I.  FUNDAMENTALS  (references  1,2.4,12.14) 

The  usual  method,  of  generating  turbulence  in  a  wind  tunnel  is  to 
place  a  grid,  or  array  of  bars  across  the  wind  tunnel  stream. 

About  15  mesh  diameters  downstream  of  the  grid,  the  turbulent  field 
is  fully  developed,  and  at  any  point  in  the  field  the  stream  velocity  can  be 
seen  to  consist  of  a  mean  component  and  a  component  fluctuating  about  this 
mean  (see  Fig.60) . 

i.e.,  U  (z,t)  =  U  (z)  +  u  (z,t)  (al) 

To  determine  the  field  of  turbulent  motion  statistically  we  require 
the  system  of  complete  joint  probability  distributions  of  the  values  of  the 
turbulence  velocity  vectors,  u(z,t),  at  any  n  points  of  space-time. 

Measurements  are  made  by  placing  hot  wire  anemometer  probes  at 
fixed  points  relative  to  the  grid  and  recording  the  variations  with  time  of 
the  linearized  output.  An  average  value  of  the  required  function  of  velocity 
is  obtained  by  taking  a  mean  value  over  a  long  period  of  time. 

If  we  assume  that  the  turbulence  signal  may  be  represented  as  a 
stationary  Gaussian  random  process  with  zero  averages  (which  it  can  if  the 
mean  value  is  subtracted  from  equation  (a  l)),  we  may  use  this  time  average 
assuming  ergodicity,  to  replace  the  probability  average  for  the  experimental 
field.  Then  all  the  relevant  statistical  information  will  be  contained  by 
the  correlation  function  (or  power  spectrum). 

As  shown  in  Appendix  B,  a  two  point  lateral  space- time  correlation 
tensor  is  required  to  evaluate  the  turbulent  field,  so  this  case  will  be 
discussed. 


Consider  two  signals  u(a,t),  u((3,t)  from  two  probes  at  O.  and  p 
(a  and  3  are  position  coordinates)  in  the  turbulent  field.  The  cross-corre¬ 
lation  for  two  such  signals  separated  in  time  by  the  lag  time  t,  is  defined 
as 

1  f ^ 

R  (d»PjT)  =  Uni  2T  I  u(a,t)  u(p,t  +  r)  dt  (a2) 

T->*°  J  _T 

where 

i  r +T 

Ruu  (a,a>T)  =  lim  2T  u(a»t)  u(a»*  +  T)  dt  (*3) 

T-»  OO  J>  '  ' 

-T 


is  called  the  autocorrelation,  and 

i  r 41 

Ruu(a>a>°)  =  lim  2T  I  u^a»t)  u(a>t)  dt 
T->  00 


-T 


is  the  variance  of  the  process. 


(a4) 


Al 


iiiii  .hi  iiii'iiimiwiniiii  mirmirw  nniinnTTf""' 


Equation  (a2)  may  also  be  expressed  as 
f  00 

Ruu  (Q,P,t)  =J  0uu  (a,p,f)  ei27rfT  df  (a5) 

-  00 

where  0uu(a,0,f),  the  cross  spectral  density  function,  is  given  by 

/T 

u(a,t)  u(P,t)  e  k7Tft  dt  (a6) 

.r 

The  function  0  (a,a,f)  is  termed  the  "power  spectrum",  and 

0Uu(O:>a:>^)  dt  represents  the  contribution  to  the  total  variance  of  the  portion 
of  the  signal  in  a  frequency  bandwidth  of  df  at  frequency  f, 

f  w  o 

i,e*’  J  \u  df  =  a  (a7) 

-  00 

The  cross  correlation  and  the  cross  spectrum  form  a  Fourier  trans¬ 
form  pair,  which  may  be  expressed  in  a  single- sided  form  as 

r  00 

Ruu  (a’P’T)  =  2  I  ^uu  (a’P’f)  gi  2jrfT  df  (a8) 

J0 

r  oo 

0uu  (a,p,f)  =  2  Ruu  (a,p,T)  e  -i2vfT  dT  (a9) 

^  o 

The  cross  correlation  function  Ruu(a,p,T)  for  u(x,t)  is  in  general 
composed  of  parts  that  are  even  and  odd  in  t 


i.e.,  Ruu  -  R^  T)  +  Rg  (^>^>T)  (aio) 


where  Rt_  and  Rg  are  as  in  Fig.6l.  The  corresponding  spectrum  function  is 


0 

uu 

(a,e,f)  = 

0j_  -  i02  (a,p,f) 

(all) 

where 

*1 

(a,&,f)  =  2 

p  00 

I  R  (a,p,r)  cos  2v  fr  dr 

J  o  uu 

(al2) 

*2 

(a,p,f)  =  2 

f  00 

/  R  (a,0,T)  sin  2 w  fT  dr 

J  o  uu 

(a!3) 

The  cross  correlation  is  a  mean  lagged  product  of  two  signals,  one 
advanced  relative  to  the  other.  Obviously,  for  a  statistically  stationary 
processes,  advancing  the  first  is  the  same  as  delaying  the  second,  so  that 
interchanging  the  order  is  the  same  as  changing  the  sign  of  t.  Accordingly 
we  have  the  reciprocity  properties 


A2 


\  (a,P»T)  =  R1  (P,a,  -t)  (a\4) 

^2  =  “^2  T)  (al5) 

The  corresponding  spectral  relations  are 

\  (<*,P,f)  =  0X  (P,  a,  f)  (al6) 

<P2  {  a,p,f)  =  02  (p,  a,  f)  (al7) 


2.  PRACTICAL  MEASUREMENT  (reference  2) 


To  exactly  determine  the  correlation  function  defined  y  (a2)  a 
perfectly  measured,  infinitely  long  record  of  random  data  is  needed.  As  this 
is  obviously  unobtainable,  and  in  fact  there  is  usually  only  a  very  limited 
amount  of  data  available  the  estimate  obtained  for  the  correlation  function 
will  differ  from  the  exact  value.  The  estimated  correlation  function 
Ruu  (ct,{3,t)  at  lag-  t,  ;>ay  be  expressed  as 


(<*,P, 


t) 


u(a,t)  u(p,  t  +t)  dt 


(al8) 


where  T  is  the  averaging  time  and  t  is  the  lag  to  be  used. 

Not  only  is  this  estimated  correlation  function  truncated  relative 
to  the  exact  function,  but  for  the  lag  times  available,  it  will  also  possess 
different  magnitudes  as  a  result  of  the  finite  averaging  time  although  due 
to  ergodicity 

!■  <RUU  (o,P,t)>  =  Ruu  (Of,p,T) 

Fortunately  it  has  been  found  that  the  estimated  correlations  may 
be  modified  by  some  even  function  of  t.  While  these  modified  estimated  corre¬ 
lation  functions  may  not  be  good  estimates  of  the  true  correlation,  their 
Fourier  transforms  can  be  very  respectable  estimates  of  "smoothed"  values  of 
the  true  spectrum  I 

The  modified  estimated  correlation  function  may  be  written 

Ruu  (a’P»  T)  =  D(T)  \n  (a>^>  T)  (a19) 


where  D(t)  may^be  regarded  as  a  window  of  variable  transmission  which  modifies 
the  values  of  Nuu  (a,p,T)  for  different  lags,  and 

D  (0)  =  1 

D  (t)  =  0  for  | t  J  >  t 
x  '  11  max 

1,  K  )  indicates  an  enseidjle  average. 
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‘■’♦rtfl’JtfU'tWfHJ*  .<*•*(  i'*,‘  "■Mil  iff'V  1 


Even  though  Ru  (o:,P,t)  was  not  defined,  for  j t  j  >  Tn,a>, ,  Ruu  (ct,r,T)  is  find  so 
has  a  definite  Fourier  transform  0^u  (a,E3,f)  which  satisfies  the  convolution 
integral 


V  ■ \  “W  fiuu  (a>e-T  1  '  'i2lr  fT 

—  00 


=  f  D(t)  f  0uu  (a,p,f)  e  l£/7f'T  df’  e-i;i77fT  dT 
J  —  00  L  sLoo 


i  .e. , 


/oo  ^ 

Q  (f-f1)  0iU  (a,p,r)  df’ 

.  m 


(aa) 


Even  though  0^u  (a,p,f')  is  not  determinate  since  R'uu(0!,P,t)  is 
not  specified  for  |t  |>t  ,  we  have 


K\n  (o,PtT)>  =  D(t)  R  (a,p,-r) 


r  00 

<<2>uu  (a,p,f)>  =  Q(f-f’)  0uu  (a,p,f’)  df 

’J-OO 


(*22) 


(*23) A. 


It  can  be  seen  that<$  (a,(3,f)^  is  an  average  over  frequency 
of  the  true  power  spectral  densityu0uu  (a,p,f'}  over  frequencies  near  f ' , 

In  general  one  would  like  to  concentrate  the  main  lobe  of  Q(f) 
near  f  -  0,  which  requires  a  flat  blocky  window.  At  the  same  time  the  lag 
window  must  be  smooth  and  gent  changing  to  reduce  side  lobe,  and  go  to 
aero  for  |t|>  Tmax,  so  a  compromise  is  necessary. 

One  simple  and  convenient  window  which  has  been  founi  to  work 
satisfactorily  is  the  "Hanning"  window,  (see  Fig. 65). 


,  COS  77T  _  t  l 

1  +  - for  t  I  <  r 

t  m 

max  -i 


=  0  for  |  t  I  >  t 

'  -  max 

For  the  case  of  data  uniformly  spaced  in  time  (as  in  this  experiment),  it  is 
most  expedient  to  first  transform  and  then  convolve  (by  equation  (a2l)) 


f  D(t)  e  _1  2irfT  dT  -  Q  (f) 

J  -00 


i.e., 


Q(f)  =  f  (I  +  k  c«s  “  )  e'l27rfT  dT 

max 

=  k  5(f)  +  £  [  0  (f+f')  +6  (f-f1) 


(a23)B. 


where  f"  =i/2t 


s  m 

0uu  (a,p,f)  =  j  jk  &(f-^)+  £  [6  (f+f'-M+  B(r.f'-A)]j0uu(a,p,M  dA 

(•24) 

Now1/2t  =  f'  =  Af  represents  the  spacing  of  the  frequency 


estimates,  i.e., 

im  (“>M)  -  tfuu  (f)  ♦  i  [  V  (f  -  Af)  t  «uu(f-Af) 


(•25) 


The  application  of  (a  25)  is  very  simple,  and  for  a  single  sided 
spectrum  (double  sided  spectra  were  considered  above)  consists  of  smoothing  in 
the  following  manner. 

J  (f  =  0)  =  .5  $  (f  »  0)  +  .5  J  (A  f) 

*uu  uu  *uu  ' 

$uu  (f)  =  .25  0UU  (f-Af)  +  .5  fjf)  +  .25  «uu(f+Af)  (a26) 

Z  (f  =  f  )  =  .5  <h  (f  )  +  .5  <£  (f  -Af) 
vuu  v  max'  '  ^uu'  max'  J  *uuv  max  ' 


Finally,  in  order  that  our  modified  spectral  estimates  exhibit  as 
much  independence  as  possible,  only  every  second  frequency  point  is  retained 
after  application  of  the  Hanning  window. 

3-  SCALE  OF  TURBULENCE  AMD  WAVE  LENGTH 

The  lengths  L  and  L  defined  by 
X  z 


Lx(z)  =  U  (z)  f 


*uu(<W) 

\u(o.“.o) 


(•27) 


■  l, 


“  Ruu 


o  R  (a,a,o) 
uu  ’ 


(•28) 


are  called  the  longitudinal  and  later  a1.  ‘ales  respectively.  They  are  convenient 


measures  of  the  linear  extent  of  tile  region  within  which  velocities  ..re  appreci¬ 
ably  correlated.  In  the  particular  case  of  shear  turbulent  flow  these  scales 
vary  with  lateral  position,  and  hence  the  notation  L(z)  where  z  is  lateral 
(vertical)  position.  ^ 

When  the  ratios  are  plotted  versus  a-P,  for  the  later' 1 

R  (a,a,o) 

scale  it  is  found  that  the  resultant  curve  depends  on  whether  (3  is  taken  upward; 
or  downwards  from  the  fixed  probe  position.  This  can  be  seer,  in  Fig. ^7.  The 
quoted  acale  values  are  averages  of  the  two  separate  scales. 

The  significant  turbulence  wave  lengths  for  excitation  of  normal 
modes  are  defined  by 


(a. 29) 


where  f  is  the  frequency  of  the  nth  mode  of  vibration  of  the  cylinder.  As  long 
as  ^  >>  b,  it  is  reasonable  to  assume  that  the  pressure  distribution  over  the 
cylinder  is  two  dimensional.  Table  2  shows  \,^2  anc*  ^3  an<*  5uKgests  that  the 
assumption  of  two-dimensionality  is  reasonable,  for  the 'rirst  modes  and  borderline 
for  the  second. 


From  the  form  of  the  narrow  band  correlation  (reference  16)  Fig.hO 
one  would  expect  that  failure  of  the  theory  would  occur  for  frequencies  where 
A  approaches  the  lateral  dimension  of  the  structure.  In  such  a  case,  you  would 
not  have  correlated  flow  over  the  cylinder,  and  the  *wo-dinensional  assumptions 
would  be  suspect. 
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APPENDIX  B;  STRUCTURE  RESPONSE 

The  following  is  a  summary  of  the  "Theory  of  the  Response  of  a 
Slender  Structure  to  a  Turbulent  Wind  with  Shear",  by  B.  Etkin.  For  the 
complete  analysis  see  reference  7- 

1.  PRELIMINARIES 

The  situation  is  shown  in  Fig.6o,  U(z)  is  the  mean  wind  profile 
(see  Fig. 2)  and  u(x,t)  is  the  fluctuating  component  (see  Fig. 3)  In  the 
direction  of  the  mean  wind  (drag  direction).  X(z)  is  the  mean  deflection  of 
the  structure  associated  with  U(z),  and  x(z,t)  is  the  vibrational  motion 
associated  with  u(z,t).  The  local  running  load  on  the  structure  is  W(z,t)  of 
which  W(z)  is  the  mean  and  w(z,t)  the  fluctuating  part.  Thus 

U  (z,t)  =  U  (z)  +  u(z,t)  (b.l) 

W  (z,t)  =  W  (z)  +  w(z,t)  (b.2) 

X  ( ?. , t )  =  X  (z)  +  x(z,t)  (b.3) 

and  the  relative  wind  is 

U’  (z,t)  =  U  (z,t)  -  y  (z,t)  (b.U) 

Assuming  the  local  load  W(z,t)  per  unit  length  is  given  by  "strip 

theory"  as 

W(z,t)  =  iL(z)  U’2(z,t)  +  i2( z)  U;2(z,t)  (b.5) 

where 

11  (z)  =  CD  (z)  b  (z)  p  /2  (b.6) 

and  „ 

12  (z)  =  k  (z)  b  (z)  p  (b.7) 

Equations  (b.5)  to  (b.?'),  as  used  subsequently,  imply  that  k  and 
C]}  are  independent  of  the  turbulence  .  This  assumption  requires  that  the 
lateral  dimension  is  small  compared  to  the  wave  length,  X,  of  the  turbulence 
defined  as 

X  =  U  /f  (see  A3)  (b.6) 

2.  FLUCTUATING  LOAD 

From  (b.5),  (b.6),  (b.7)  and  (b.8)  we  have 
W(z,t)  =  ^(zHU2-  2  Ux  +  x2  ]  +  i2(z)  [U-  x*]  (b.9) 

Writing  the  fluctuating  load  as 

w  (z,t)  =  W(z,t)  -  W  (z)  (b.10) 
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where 


gn(z)  =  2  U(z)  i^z)  Fn(z)  (b.19) 

hn(z)  =  i2(z)  Fr(z)  (b.20) 

Figure  62  shows  the  manner  in  which  the  individual  modes  are  excited  and 
contribute  to  the  total  deflection  x. 

1* .  SPECTRAL  ANALYSIS 


The  spectral  density  of  x  is  given  by  the  fundamental  equation 

0  (z,f)  =  £  Z  F  (z)F  (z)  G*  (if )G  (if)0„_(f)  (b.2l) 

xx  ’  m  n  m'  n  mv  '  n 

is  the  cross  spectral  density  ofi(t)  and^7n(t),  while  Gn(if)  is  the 
transfer  function  relating  x(z,t)  to “^(t). 

From  the  equations  of  motion 

Gn(if)  “  [In(  M \ 2-  +  8  It r2^  ff  )  ]-1  (b.22) 

n 

The  spectrum  function  is  the  basic  information  needed 

to  assess  the  behaviour  of  the  structure. 


To  obtain  the  general  spectrum  functions  for  the  driving  forces 
in  terras  of  the  turbulence  velocity  spectra,  we  proceed  via  the  correlation 
function. 


<T) 


-i27rft 

e 


dT 


(b.23) 


where 

V<T)  “<3m't)3n  (t  T» 


and^,  y represents  an  ensemble  average  which  is  equivalent  to  a  time  average 
for  this  case. 


On  forming  the  appropriate  products  of  then's  (given  by  equation 
(b.l8))  as  shown  by  equation  (b.2^)  and  noting; 

(i)  that  the  order  of  integration  and  averaging  can  be 
interchanged 

(ii)  that  the  mean  products  in  these  integrals  are  them¬ 
selves  cross  correlations 


(iii)  that 


0 .  ,  =  (J  0 

uu  uu 


0.  =  -ito  0 

uu  uu 
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we  obtain  the  cross  correlation  of  *3  and 3  >  which  in  turn  can  be  Fourier 

fll  n 

transformed  to  give  the  required  spectrum,  The  resultant  generalized  force 
spectrum  is 

%^(f)  =  //  0uu(a,P,f){gn(a)g m(P)  +  **ir^ 
o 


+  i  2mf[gn(a)hm(p)-  gm(P)hn(Q) ]|  da  d£  (b.25) 


Obviously  the  only  spectrum  function  needed  to  evaluate  all  the  terras  is 
0uu(a,p,f),  the  velocity  cross  spectrum. 


APPENDIX  C:  SYSTEM  PARAMETERS 
1 ■  GENERAL 

m  =  mass/unit  length  =  2.080  x  10  ^  slugs/ft. 
k  -  2-D  added  mass  coefficient  =  tt/4 
CD  -  1.2 

l  =  total  cylinder  length  =29*5  inches 
1'  =  exposed  cylinder  length  =  27.0  inches 
b  =  cylinder  diameter  =  .50  inches 
p  =  2.38  x  10~"  slugs/ft. 

p 

EI  -  bending  stiffness  =  518O.O  lb. inch 
m  =  m+kb  p  =  2.083  x  10  J  slugs/ft. 


2.  MODE  FREQUENCIES  AND  SHAPES 


vibration  is 


The  differential  equation  for  a  free  slender  beam  undergoing  small 


(EI  x")  +  mx  =  0 


(c.l) 


where  x  =  x(z,t)  is  the  lateral  deflection.  By  applying  boundary  conditions 
applicable  to  a  cantilever,  it  can  be  shown  (reference  8)  that  the  natural 
mode  frequencies  are  given  by  (see  table  l) 

a  2  r  l 

f  _  n  £.  S  2  in  o\ 

f"  '  T  J  L“  J  '  } 


where  =  .597,  =  1.49,  =  5/2>  a4  =  7/2  . 

The  natural  mode  shapes  are  given  by  (see  Fig. 63), 

„  t  \  sin  pi-  sinh  pi]/.,  ,  \ 

F  (z)  =  - n - - - t~-r-  (sinh  pz  -  sin  pz) 

n  cosn  p^>  +  cos  p^  ' 


+  (cosh  pz  -  cos  pz) 


(c.3) 


where 


■  1  2  tL  ,  q 

4  7r  *n  m 

[ - eT - _ 


The  second  derivative  of  the  mode  shape  function  evaluated  at  z  =  0  is 


Cl 


. .  mimirtmiii>wiBWilllMW|iiwoll*lwmwWW|W!,l,»fflWHd?i^ "W*- "|41 


i 


dz  |z=0 


=  F  "  (o)  =  2 
n 


2tt£  m1  i  2 
n 


i.  GENERALIZED  INERTIA 


(c.M 


The  expression  for  the  generalized  inertia  is  given  by  equation 
(b.17)  (see  table  1). 


i.e.. 


r£ 

[  -  F 

n  Jn  n 


(z)  m’  (z)  dz 


ft  2 

I  =  m’  F  (z)  d: 
n  f  n 


( c .  5) 


for  m'  constant. 

4.  DAMPING 

(a)  Structural: 

The  structural  damping  was  measured  for  the  first  3  modes  by 
recording  the  decay  with  time  of  the  signal  from  the  cylinder  strain  gauges 
on  visicorder  tape  after  the  signal  was  passed  through  a  narrow  filter  centered 
about  the  natural  frequency  in  question.  An  exponential  decay  scheme  repre¬ 
sentative  of  a  second  order  system  was  assumed,  i.e., 


v  -t  2  itf  t 

i(t)  =  eft  e  isn  n 


(c.6) 


where  eQ  was  the  initial  signal  and  £Sn  is  the  structural  damping  ratio. 
From  (c.6) 


in  (e/eQ)  =  -  £  2  7/f  t 

'  '  °  Ss  n 

n 


thus  (see  table  1) 


tL  =  (  -l/2vfnt)  in  (e/eD) 


(c.7) 


The  variation  of  damping  with  initial  displacement  was  also  investigated  and 
the  results  are  shown  in  Fig. 64. 

(b)  Aerodynamic: 

The  aerodynamic  damping  coefficient  per  unit  length  is  given  by 


cjz)  »  2  ix  U  (z)  F n(z)‘ 


(c.8) 


Thus  the  damping  ratio  for  the  total  model  length  may  be  expressed  as  (see 
table  1) 

Si  r  S' 

?a  =  iV  I  »„(*>'"  0  (»)  d.  (c.9) 

n  n  n  ^  o 


(c)  Total: 


The  total  damping  is  given  by  (see  table  l) 
^T  =  ^a  +  ^s 


n  n  n 


5.  TRANSFER  FUNCTION 


The  transfer  function  for  the  second  order  system  is  ( from (b. 22) 


(2  77fn)2  -  (277f)2  +  i  8/f  f  CT  \ 

n->  J 


and  so 


|  Gn(if)|2  =  |ln2  [Wf^-  ^2)  +  S^ttVCj  2fn2  | 


Equation  (c.ll)  is  shown  in  Fig. 66. 


APPENDIX  D:  ERRATA  FOR  DEVIOUS  CAMPBELL  *  ETKIN  REPORT 


1.  Appendix  A:  Equation  (21),  p.A6  should  read 
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1  r 00 
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2.  Figure  10:  The  ordinate  designation  should  read 
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TAPLE  3 


Mean  velocity 


Fin. 2  -  MEAN  VELOCITY  PROFILE 
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Fig. 3  -  TURBULENCE  INTENSITY  SURVEY 


Fig. 4  -  GENERAL  VIEW  OF  MODEL 


Fig. 5  -  STRAIN  GAUGES  MOUNTED  ON  MODEL 
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Fig. 6  -  MODEL  MOUNTING  STRUCTURE 


MOUNT  INSTALLATION  DETAIL 


Fig. 9  -  WIND  TUNNEL  AERODYNAMIC  OUTLINE 


g(f/f 


ANALYSER  FILTER  BAND-SHAPE 
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Strain  Gauge  Signal  Analysis 


Hot  Wire  Linearizing  Circuit 
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Fig. 14  -  HOT  WIRE  CALIBRATION 


Fig. 15  -  HOT  WIRE  PROBES 


Figs. 16-19  -  SPACE-TIME  CORRELATION  COEFFICIENTS 


Figs. 20-22  -  SPACE-TIME  CORRELATION  COEFFICIENTS 


Figs. 23-26  -  SPACE-TIME  CORRELATION  COEFFICIENTS 
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Fig. 27  -  REYNOLD'S  NUMBER  DISTRIBUTION 
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Figs. 28-32  -  TURBULENCE  CROSS  SPECTRA 


Fig. 33  -  TYPICAL  TURBULENCE  POWER  SPECTRUM 


Figs. 34-38  -  TURBULENCE  CROSS  SPECTRA 
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Fig. 39  -  GENERALIZED  INPUT  FORCES 


.40  -  NARROW  BANO  CORRELATION  FUNCTION 


Figs. 41-46  -  NORMALIZED  CROSS-CORRELAT JON  FUNCTIONS 


Figs. 47-51  -  NORMALIZED  CROSS-CORRELATION  FUNCTIONS 


Fig. 52  -  LATERAL  INTEGRAL  SCALE 


Fig. 53  -  LONGITUDINAL  INTEGRAL  SCALE 
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Fig. 54  -  STROUKAL  VORTEX  SHEDDING  FREQUENCY 


Fig. 55  -  EXPERIMENTAL  LIFT  BENDING  MOMENT  SPECTRAL  RESPONSE 
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Fig. 56  -  COMPARISON  OF  FILTERED  THEORETICAL  AND  EXPERIMENTAL  SPECTRA 
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Fig. 58  -  SPECTRUM  COMPARISON 
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Fig . 59  -  NARROW  BAND  VARIANCES 
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Fig. 61  -  SYMMETRIC  ANO  ANTISYMMETRIC  COMPONENTS  OF  THE  CORRELATION  FUNCTION 


Fig. 62  -  SUPERPOSITION  OF  MODEL  RESPONSES 


Fig. 63  -  CANTILEVER  CYLINOER  MODE  SHAPE  FUNCTIONS 
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Fig. 64  -  STRUCTURAL  DAMPING 
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Fig. 67  -  C  Y  L.  t  N  fit  R  STRAIN  GAUGE  OUTPUTS 


